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Abstract 



We derive rigorous upper and lower bounds for the ground state entropy of the g-state 
Potts antiferromagnet on the honeycomb and triangular lattices. These bounds are quite 
restrictive, especially for large q. 
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Nonzero ground state disorder and associated entropy, So 7^ 0, is an important subject 
in statistical mechanics; a physical realization is provided by ice, for which So = 0.82 ± 0.05 
cal/(K-mole), i.e., So/ks = 0.41 ± 0.03 [l], |J. Ground state (g.s.) entropy may or may not 
be associated with frustration. An early example with frustration is the Ising (equivalently, 
q = 2 Potts) antiferromagnet on the triangular lattice |§. However, g.s. entropy is also 
exhibited in the simpler context of models without frustration, such as the g-state Potts 
antiferromagnet (AF) f4|-|| on the square (sq) and honeycomb (he) lattices for (integral) 
q > 3 and on the triangular (tri) lattice for q > 4. Of these three 2D lattices, Sq has been 
calculated exactly for the triangular case @, but, aside from the single value So(sq,q = 
3)/fcs = (3/2) ln(4/3) ||, not for the square or honeycomb lattices. Therefore, it is valuable 
to have rigorous upper and lower bounds on this quantity. Using a "coloring matrix" method, 
Biggs derived such bounds for the square lattice ||. Here we shall extend his method to 
derive analogous bounds for the honeycomb lattice and compare the results with our recent 
Monte Carlo measurements [[H], [TTJ] and with large-g series [12]j. We also derive such bounds 
for the triangular lattice; the interest in this case is that the bounds can be compared with 
the exact result of Baxter j7| . 

We make use of the fact that the partition function at T = 0, Z(A,q,K = — oo), for 
the g-state zero-field Potts AF on a lattice A (where K = (3 J, (3 = l/(fceT), and J < 
denotes the spin-spin coupling) is equal to the chromatic polynomial P(A, q). Here, P(G, q) 
is the number of ways of coloring the graph G with q colors such that no adjacent vertices 
(sites) have the same color [13|]. Define the reduced, per site free energy for the Potts AF 
in the thermodynamic limit as f(A,q,K) = limjv^oo AT -1 In Z(A, q, K). From the general 
relation between the entropy S, the internal energy U, and the reduced free energy, S = 
/3U + f (henceforth, we use units such that k,B = 1), together with the property that 
\\m K ^_ QO (3U ((5) = 0, as is true of the g-state Potts AF models considered here, it follows 
that So (A, q) = f(A, q,K = — oo) = In W(A, q), where W(A, q) is the asymptotic limit 

W(A,q)= lim P(A,q) 1 / N (1) 

N^oo 

Given this connection, we shall express our bounds on the g.s. entropy S (A, q) in terms of 
the equivalent function W(A, q). As we have discussed earlier [|llj, the formal eq. ([!]) is not, 
in general, adequate to define W(A, q) because of a noncommutativity of limits 

lim lim P(A,q) 1/N + lim lim P(A,q) 1/N (2) 

at certain special points q s . We denote the definitions based on the first and second orders of 
limits in (H) as W(A,q)c nq and W(A,q)r> qn , respectively. This noncommutativity can occur 
for q < q c (A), where q c (A) denotes the maximal (finite) real value of q where W(A,q) is 
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nonanalytic Since q c (hc) = (3 + Vo)/2 = 2.618... and q c {tri) = 4 (and q c (sq) = 3) [|TT 



it follows that for the ranges of q considered here, viz., (real) q > 3 for A = he, sq and q > 4 
for A = tri, one does not encounter the noncommutativity (0). 

To proceed, we consider a sequence of (regular) m x n 2D lattices of type A, with 
periodic boundary conditions (PBC's) in both directions, and m and n even to maintain the 
bipartite property for finite square and honeycomb lattices and thereby avoid frustration. 
For A = sq, Biggs introduced the notion of a coloring matrix T, somewhat analogous to 
the transfer matrix for statistical mechanical spin models. The construction of T begins 
by considering an n- vertex circuit C n along a column of A, i.e., a ring around the toroidal 
lattice, given the PBC's. The number of allowed g-colorings of this circuit is P(C n , q) = 
(q— l)[(q — l) n_1 + (—1)™]. Now focus on two adjacent columnar circuits, C n and C' n . Define 
compatible g-colorings of these adjacent circuits as colorings such that no two horizontally 
adjacent vertices on the circuits have the same color. One can then associate with this 
pair of adjacent columnar circuits an M x M dimensional symmetric matrix T, where M = 
P(C n , q) = P{C' n , q) with entries T Cn ,c n = T c < n ,c n = 1 or if the g-colorings of C n and C' n 
are or are not compatible, respectively. Then for fixed m, n, P(A mxn ,q) = Tr(T m ). For a 



given n, since T is a nonnegative matrix, one can apply the Perron- Frobenius theorem ]14 
to conclude that T has a real positive maximal eigenvalue \ ma x,n(q)- Hence, for fixed n, 

Jim Tr^™) 1 /^ - \]i: x (3) 

so that 

W(K,q)=\m L X 1 J i : x (4) 

Denote the column sum Kj(T) = Y^LiTij (equal to the row sum Pj(T) = YnLiTji since 
T T = T) and S(T) = E^=i% note that S(T)/Af is the average row (column) sum. 
Combining the bounds for a general nonnegative J\f x M matrix A, flT4 



min{7 j (A)} < \ max {A) < max{7 j (A)} , for ^ = Kj or pj (5) 

with the (k = 1 case of the ) more restrictive lower bound applicable to a symmetric non- 
negative matrix 15|| , 

-S(T k )i 1 / k 



we have 



S(T) 



< X max , for k = 1,2, ... (6) 
< ^max{T) < maxl^T)} (7) 



N 

For the honeycomb lattice, we find that the analogue of the circuit C n on the square 
lattice is the set of vertical dimers shown in Fig. 1(a), which we denote as p. With m and 
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n even to maintain a bipartite lattice, there are n/2 dimers in p, and the total number of g- 
colorings of these dimers is A^ C) „ = [q(q — l)] n ^ 2 . We next associate the T matrix T(hc n ) with 
two adjacent sets of dimers p and p' (see Fig. 1(a)); T{hc n ) is thus a Nh c , n x ^V/ic.n matrix. 
Two g-colorings of the dimer sets p and p' are compatible if and only if the horizontally 
adjacent vertices have different colors, and T p y = 1 or 0, respectively, if these colorings are 
compatible or incompatible. We observe that S(T) = -P(C2„,g). Therefore, 

Spike*)) (q-l^q-iy^ + l 



N hc , n [q(q ~ l)]"/ 2 

To calculate the maximal column sum, we consider two neighboring sets of dimers p 
and p', with n sites each labeled by i = 0, 1, ...,n — 1 (see Fig. 1(a)). Let the sites of set 
p be colored in such a way that sites on the same dimer have different colors (choosing 
one such configuration of colors corresponds to fixing one column in the color matrix T). 
Let Xj denote the number of g-colorings of sites i = to i = j of p', such that a site i 
in set p has a different color from a site i on p'. If j is odd, the coloring of the j'th site 
in p' is only constrained to be different from the coloring of the adjacent j'th site in p, so 
Xj = (q — l)Xj-i. |TJ| The color assigned to an even-j site in p' must be different from the 
color of (i) the other member of the dimer in p' and (ii) the adjacent j'th site in p; hence, 
Xj = {q — 2)Xj-_i + Yj-i, where Yj-\ denotes the number of colorings for which site j of p 
has the same color as site j — 1 of p'. Note that Yj-\ is a subset of Xj-%, i-e. Y,-i < 
Thus 



X 3 <{q- 2)X j . 1 



X 



j-2 



Xj = {q- 1)X 



Using eq.(|l0|) in eq.(||) and setting j 



X21 < 



(g-l)(g-2) + l 



for even j, 
for odd j, 
2£, we have 

X 2( -2 , for 



which yields X 2 ^ < (q — l)(q — 2) + 1 X , where X = g 



X n _x = (g - 1)X„_ 2 <(q- lY[(q - l)(g - 2) + 1 

Because max{/«j(T(/ic n )} < X n _ 1; we obtain 

max{fi:(T(/ic„))} < (g - l) 2 [g 2 - 3g + 3] 

Hence, using (|3|) and (ffl), we derive the bounds 
(g - I) 3 / 2 



2 < j < n - 2 

1 < j < n - 1 

2 < 2£ < n - 2 , 

— 1. It follows that 

(n-2)/2 



(n-2)/2 
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1/2 



< H/(/tc,g) < (g 2 -3g + 3) 1/2 for g>3 



(9) 
(10) 

(11) 

(12) 
(13) 

(14) 
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Figure 1: (a) Honeycomb and (b) triangular lattices. See text for discussion. 



The bounds are also seen to apply for the case q = 2 if one uses as a definition W(hc, 2) = 
W(hc, 2) Dnq = 1 given by the first order of limits in eq. (H). 

Similarly, for the triangular lattice we define the color matrix T(tri n ) by considering the 
compatibility of g-colorings of two neighboring n-vertex circuits C n and C' n . An example of 
such adjacent circuits is shown by the darker lines in Fig. 1(b). Here T(tri n ) is a N tr i, n x N tr i^ n 
matrix, where Ntn, n = P(C n , q) = P(C' n , q). For the triangular lattice with periodic or open 
boundary conditions in the vertical direction, S(T) is equal, respectively, to the numbers 
P(cctri n , q) and P(octri n , q) of g-colorings of a cyclic or open chain of triangles with 2n ver- 
tices. In the n — > oo limit of interest here, lirn n _»oo P(cctri n , q) 1 ^ = lim n _ +QO P(octri n , q) 1 ^, 
so it does not matter which type of chain we use. An elementary calculation yields 



P(octri n , q) = q(q - l)(q - 2) 



2n-2 



(15) 



so 



lim 

n— >oo 



S(T{tri n )) 



l/n 



(g ~ 2) ; 
q-1 



(16) 



Ntri,n 

To calculate max{Kj(T(tri n ))}, we derive, as before, an upper bound for X n -\. Each 
vertex of C' n is connected to two vertices of C n , hence each of the Xj-\ colorings of the sites 
i = to % = j — 1 of the circuit C' n can be extended in at least q — 3 ways to the site i — j. 
Thus, for the triangular lattice, the equivalent of eqs. 



and (1C) is 



Xj <(q- 3)X j _ 1 + X 



i-2> 



2 < j <n- 1, 



(17) 



with Xo = q — 2 and X\ < (q — 3)(q — 2) + 1 [17|. This recursion relation is of the same 
form as the one obtained previously for the square lattice, and can thus be solved by the 
same method. We thus find for the triangular lattice the inequality 

(<7-2) 2 



q 



! < W(tri, q) < ]- [q - 3 + {q 2 - Qq + 13) 1/2 



for q > q c (tri) = 4 



(For q = 3, the ground state entropy is zero, i.e., W(tri, 3)D m = 1-) 

Denote the lower and upper bounds for these lattices as W(A, q)i and W(A, q) u , respec- 
tively. We observe that as q increases, these bounds rapidly approach each other, and hence 
restrict the exact values W(A, q) very accurately. This can be seen as a consequence of the 
fact that, aside from the obvious prefactor q, W(A, q)e and W(A, q) u are the same up to 
0{q- 2 ): 



q- 1 W(hc,q) i 
q- l W(hc,q) u 



:1 



8 



q 



l 

16 



q- 3 + 0(q-*) 



l_ ! g -i + 3 - 2 + -g- 3 + 0(g" 4 ) 
2 H 8 y 16 y KH ' 



(19) 
(20) 
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Q 


W(hc,q) e 


W(hc,q) u 


W(sq,q) e 


W(sq,q) u 


W(tri,q)i 


W(tri,q) u 


W(hc,q) MC 


W(hc,q) MC 


W(sq,q) 


W(sq,q) 


W{tri,q) 


W(tri,q) 


3 


0.98390(60) 


1.04358(65) 


0.97425(55) 


1.05091(60) 


— 


— 


4 


0.99781(60) 


1.01612(60) 


0.99844(65) 


1.03305(65) 


0.91262 


1.107485 


5 


0.99948(55) 


1.00726(55) 


0.99970(60) 


1.01593(60) 


0.99377 


1.06630 


6 


0.99978(65) 


1.00377(65) 


0.99992(60) 


1.00851(60) 


0.99879 


1.03087 


7 


0.99988(65) 


1.00220(65) 


0.99996(60) 


1.00498(60) 


0.99963 


1.01628 


8 


0.99999(60) 


1.00145(60) 


0.99996(65) 


1.00312(65) 


0.99986 


1.00953 


9 


1.00001(60) 


1.00099(60) 


0.99995(65) 


1.00206(65) 


0.99994 


1.00602 


10 


0.99994(60) 


1.00063(60) 


0.99986(60) 


1.00134(60) 


0.99997 


1.00404 



Table 1: Comparisons of lower and upper bounds with Monte Carlo measurements of 
W(hc, q) and exact values of W(tri,q), respectively, for g c (A) < q < 10. An analogous 
comparison is included for A = sq. 



and similarly with q~ l W(tri, q) b , b — £,u. (This was also true for the A = sq bounds ||). 

In Table 1, we compare the bounds (|I~4D for A = he with our recent Monte Carlo measure- 
ments of W(hc, q) |T0, II]. We also compare our bounds ( Jl8| ) for A = tri with the exactly 



known results of Baxter J7). For reference, Table 1 includes a similar comparison of the 



A = sq bound with the known q = 3 value || and Monte Carlo measurements |18|, [H], |Tl]] for 
q > 3. We see that as q increases past q ~ 4, the upper and lower bounds bracket the actual 
respective values quite closely, and that the latter values lie closer to the lower bounds. 

To understand why the actual values of W(A, q) lie closer to the respective lower bounds, 
we compare the large-g series with the expansions of these lower bounds. For a lattice A 
with coordination number (, the large q series can be written in the form 

W(A,q)=q(l-q- l )^ 2 W(A,y) (21) 

where W(A,y) = 1 + J2™=i w n y n with y = Defining the analogous functions W(A,y) b 
via 

W(A, q) b = q(l - q- l f/ 2 W(A, y) b , b = £, u , (22) 

we obtain W(hc,y)e = 1, which agrees to the first five terms, i.e., to order 0(y 4 ), with the 
series W{hc, y) = 1 + y 5 + 2y n + 0(y 12 ), while W{hc, y) u = l + y 3 /2 + 0(y 6 ). We also 
calculate W(tri,y)i = (1 — y 2 ) 2 , which agrees to the first five terms, i.e. to 0(y 4 ), with the 
series expansion of the exact Baxter result, W(tri,y) = 1 — 2y 2 + y A + y 5 + 0(y 6 ), while 
W(tri,q) u = 1 — 2y 2 + 2y 3 + 0(y 4 ). Finally, W(sq,y)e = 1 + y 3 , which agrees to the first 
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seven terms, i.e, to 0(y 6 ), with the series |L2|, [19| W(sq, y) = 1 + y 3 + y 7 + 0(y s ), while 
W(sq,y) u = l + 2y 3 + 0(y 4 ). 

In summary, we have derived rigorous upper and lower bounds for the (exponent of the) 
ground state entropy of the Potts antiferromagnet on the honeycomb and triangular lattices 
and have shown that these are very restrictive for large q. Since nonzero ground state 
entropy sheds light on some of the most fundamental properties of statistical mechanics, it 
is of interest to derive similar bounds for other lattices; work on this is in progress. 

This research was supported in part by the NSF grant PHY-93-09888. 
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